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Abstract. We use the elliptic reconstruction technique in combination with 
a duality approach to prove a posteriori error estimates for fully discrete back- 
ward Euler scheme for linear parabolic equations. As an application, we com- 
bine our result with the residual based estimators from the a posteriori esti- 
mation for elliptic problems to derive space-error estimators and thus a fully 
practical version of the estimators bounding the error in the Lcc (0, T; L2 (H)) 
norm. These estimators, which are of optimal order, extend those introduced 
by Eriksson and Johnson IEJ91I by taking into account the error induced by 
the mesh changes and allowing for a more flexible use of the elliptic estima- 
tors. For comparison with previous results, an application of our abstract 
results using residual estimators is provided. 



1. Introduction 

A posteriori error estimators and their use to derive adaptive mesh refinement al- 
gorithms to solve time-dependent problems constitute the object of intense research. 
The problem is appealing for the theorician as a test ground for novel analytical 
techniques as well as for the practitioners which are interested in minimizing the 
ammount of computational time in order to obtain a satisfactory accuracy in the 
computer simulations of time-dependent PDE's. Both the theorical and practical 
aspects of a posteriori-based adaptive numerical methods for evolution partial dif- 
ferential equations has benefitted immensely from the surge in the production of 
dedicated papers in the last 15 years. 

In a recent article jLM06j we have introduced the elliptic reconstruction for fully 
discrete explicit Euler schemes and we have applied it using the energy techniques to 
derive optimal-order a posteriori residual-based error estimators. In particular we 
were able to establish estimates for the error in the Loo(0, T; L2(r2)) without the use 
of duality methods. Previous work for the spatially-discrete and time-continuous 
scheme was introduced by Makridakis & Nochetto |MN03j . 

Our chief goal in this paper is to explore the possiblity of using the elliptic re- 
construction techinque in conjunction with the duality technique as introduced by 
Eriksson & Johnson |EJ91j . The duality technique provides an important alterna- 
tive to energy techniques and is widely used for the derivation of a priori and a 
posteriori error estimates both for elliptic and parabolic problems. Since being first 
considered by Eriksson and Johnson in 1991 [EJ91| it has been developed in many 
different directions, including its use in implicit and goal oriented a posteriori error 
estimates. 
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The elliptic reconstruction has been used in combination with energy estimates, 
in which one mimicks the energy estimates for the parabolic equation in order to de- 
rive error estimates from a PDE where the error, or part thereof, is the "unknown" . 
Our chief objective in this paper is to exhibit the flexibility of the elliptic recon- 
struction technique by showing that it can be completely decoupled from energy 
considerations. This is not obvious, indeed, in many works a posteriori analysis, the 
elliptic part is entangled with the parabolic part and there is not a clear cut differ- 
ence between elliptic and parabolic effects. As noted in recent work on a posteriori 
analysis for time-dependent problems [AMNOS) IBBMOSj IBV04|, e.g.], understand- 
ing the splitting between the elliptic, stationary, and parabolic, time-dependent, 
errors, as well as the part of the error where these effects are coupled, is important 
in designing adaptive methods. 

An important by-product of our approach is that the mesh-change in time is 
considered as part of the proofs of our theorems. Indeed, unlike former derivations 
a posteriori error estimates via duality |EJ91j . we do not impose on the mesh any 
assumption that are susceptible of violation in a practical implementation of the 
scheme, such as the no-refinement assumptions. We stress that once the elliptic 
reconstruction paraphernalia is introduced the analysis becomes quite clear and 
straightforward and the heart of the matter does not exceed more than 3 pages, in 
contrast with the more involved approach of Eriksson & Johnson. 

To justify completely our effort from a practical side, we give an application 
of our theory with a concrete example where residual-based estimators are used 
[AOOOl [LN03| . We emphasize, however, that our results are not limited to the use 
of residual-based estimators and that other estimators which work for the L2 norms 
in elliptic problems could be used. 

1.1. Basic set-up. We start by introducing the exact PDE whose discretization 
is the object of this paper. Let be a bounded domain of the Euclidean space R'^, 
d G Z+ and T G K"*". We assume throughout the paper that is a polygonal 
convex domain, noticing that all the results can be extended to certain non-convex 
domains, like domains with reentrant corners in d = 2. 

Given a Lebesgue measurable set D C M'', we use the following notation 

(1.1) (0,1^)^ := / cj){x)^P{x)d^l{x), 

(1.2) ||0|Iz,:=II0!Il.(D):=(0.0>d'. 

(1.3) := II 0*^011^ , for fc e Z+ (with DV := V0, etc.), 

/ \i/2 

(1-4) U\kD--=(Ml + T.\^\iD) ,forfteZ+, 

where dfi{x) is either the Lebesgue measure element da;, if D is has positive such 
measure, or the (d — l)-dimensional (Hausdorff) measure ds(a;), if D has zero 
Lebesgue measure. In many instances, in order to compress notation and when there 
is no danger of engendering confusion, we may drop altogether the "differential" 
symbol from integrals. This convention applies also to integrals in time. 

We will use the standard function spaces L2iD), U'^iD), UoiD) and denote by 
H^^(£>) the dual space of IlJ(£') with the corresponding pairing written as (• | •)^. 
We omit the subscript D whenever D = U. We denote the Poincare constant 
relative to by 6*2,1 and take |-|^ to be the norm of HJ(0). We use the usual 
duality identification 

(1.5) Hj(f7)cL2(J7)-L2(J7)'cH-i(f^) 
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and the dual norm 

(1.6) ||^,|L,:= sup Mi)f= sup 

Let a be the elliptic bilinear form defined on HQ(il) by 

(1.7) a {v, x/j) ■= {AVv, W) , Vw, V e Hj(f]) 

where "V" denotes the spatial gradient and the matrix A £ Loo(f^)'^^'* is such that 

(1.8) a(V;,0) </3|^U0li, V0,V^gH1(17), 

(1.9) ai^,^)>a\<f>\l, V0eHi(f}), 
with a, P Cz M+. We also use the energy norm \ -\^ defined as 

(1.10) 101, :=a(0,0)l/^ V0eHi(f^). 

It is equivalent to the norm on the space HQ(f7), in view of p.Sp and p.9p . In 
particular, we will often use the following inequality 

(1.11) |0|i<«-i/2|^|^,V(/)GHi(f}). 

Let u e Loo(0,T;Hj(Sl)), with dtu e L2(0, T; H"\f7)), be the unique solution of 
the linear parabolic problem 

and u(0) = g, 

where / € L2{n x (0, T)) and g e Hj(r2). 

Whenever not stated explicitly, we assume that the data f,g,A and the solution 
u of the above problem are sufficiently regular for our purposes. 

In order to discretize the time variable in ()1.12p . we introduce the partition 
= to < ^1 < ••• < tN = T of [0,T]. Let /„ := (t„_i,t„] and we denote by 
Tn tn — tn-1 the time steps. We will consistenly use the following "superscript 
convention": whenever a function depends on time, e.g. f{x,t), and the time is 
fixed to be t = i„, n S [0 : iV] we denote it by /"(jc). Moreover, we often drop the 
space dependence explictly, e.g, we write f{t) and /" in reference to the previous 
sentence. 

We use a standard FEM to discretize the space variable. Let {^n)ne[0:N] be a 
family of conforming triangulations of the domain O |BS94[ ICia78| . These trian- 
gulations are allowed to change arbitrarily from a timestep to the next, as long 
as they maintain some very mild compatibility requirements. Our use of the term 
"compatibility" is precisely defined in |LM06| : it is an extremely mild requirement 
which is easily implemented in practice. 

For each given a triangulation we denote by /i„ its meshsize function defined 

as 

(1.13) hn{x) ^ diam(if), where K £ 3n and x G K, 

for all a; e il. We also denote by the set of internal sides of 5^, these are edges 
in c? = 2 — or faces in d = 3 — that are contained in the interior of ft; the interior 
mesh of edges E„ is then defined as the union of all internal sides Ue£^„^- We 
associate with these triangulations the finite element spaces: 

(1.14) V" := e Hj(r!) : VA^ e ^„ : G P^} 

where is the space of polynomials in d variables of degree at most £ G Z+. 
Given two successive compatible triangulations ■%i-i and we define /i„ := 
max(/i„, hn-i) [LM06|, Appendix]. We will also use the sets S„ := E„ n and 
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1.2. Definition (fully discrete scheme) We consider the following fully discrete 

scheme of problem ()1.12p associated with the finite element spaces V": 

(1.15) 

C/° := l°u{0), and 



V$„ e V",for ne\l:N] 



Here the operator /° is some suitable interpolation or projection operator from 
Hj(f2), or L2(f2), onto V", and /" equals either the value of / at t„, /" := f{-,t 



backward (or implicit) Euler-Galerkin finite clement scheme jTho97| . 



or its time-average on /„, " ^ /(•,t)dt/(t„ — t„-i). This schme is the standard 



In the sequel we shall use a continuous piccewise linear extension in time of 
the sequence (f/") which we denote by U(t) for t G [0, T] (see ^2.41 for the precise 
definition). 

1.3. A posteriori estimates and reconstruction operators. Suppose we as- 
sociate with U an auxiliary function w : [0, T] Hj(ri), in such a way that the 
total error 

(1.16) e:=U-u 
can be decomposed as follows 

(1.17) e = p-t 

(1.18) e := ti) — [/, p := uj — u. 

The new auxiliary function uj is reconstructed from the given approximation U. 
The success of this splitting for the estimates to follow rests in the following prop- 
erties: 

1. The error e is easily controlled by a posteriori quantities of optimal order. 

2. The error p satisfies a modification of the original PDE whose right-hand 
side depends on e and U . This right-hand side can be bounded a posteriori in an 
optimal way. 

Therefore in order to successfully apply this idea we must select a suitable recon- 
structed function lu. In our case, this choice is dictated by the elliptic operator at 
hand; the precise definition is given in §2.21 In addition the effect of mesh mod- 
ification will reflect in the right-hand side of the equation for p. As a result of 
our choice for u we are able to derive optimal order estimators for the error in 
Loo(0,r;L2(f7)), as wefl as in Loo(0, T; HJ(0)) and H^(0, T; L2(f^)). In addition, 
our choosing uj as the elliptic reconstruction will have the efi^ect of separating the 
spatial approximation error from the time approximation as much as possible. We 
show that the spatial approximation is embodied in e which will be referred to as 
the elliptic reconstruction error whereas the time approximation error information 
is conveyed by p, a fact that motivates the name main parabolic error for this 
term. This "splitting" of the error is already apparent in the spatially discrete case 
|MNn3| . 

With the above notation, we prove in the sequel that p satisfies the following 
variational equation. 

1.4. Lemma (main parabolic error equation) For each n £ [1 : N], and for each 

{dtp I 0) -f a (p, 0) = {dte, 0) + a (c^ - c^", 0) 
^^■^^^ + (P^r - /, 0) + r-' {P^U--' - on /„. 

Here Pq denotes the -projection into V". 
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1.5. Outline. The rest of this article is organized as follows: in fJSJwe introduce the 
main tools related to the elliptic resconstruction, in §3 we outline the strategy of the 
proof in section §4 we present an application using the residual-based estimators. 

2. The discrete scheme and the elliptic reconstruction 

We introduce the numerical schemes that we study, some basic tools including 
the definition of the elliptic reconstruction. To keep this section brief, and due to 
a substantial overlap, we refer the reader to our previous paper for the details []. 

Although the definitions some in parts of this section are independent of the 
time discretization and could be applied to any finite element space, we still use 
the space V" defined in the introduction. 

2.1. Definition (representation of the elliptic operator, discrete elliptic operator, 
projections) Suppose a function v e V", the bihnear form can be then represented 
as 

(2.1) a{v,c^)= {-^w{AVv)A)k+ E {J[v]A) e ^^H^)^ 

where J[v\ is the spatial jump of the field AVw across an element side E G S^n 
defined as 

(2.2) J[v\\e {x) = {AVvlj^ {x) := lim {AVv{x + £Ue) - AVi'(a; - eve)) ■ ue 

where i^^; is a choice, which does not influence this definition, between the two 
possible normal vectors to E at the point x. 

Since we use the representation (|2.ip quite often, we introduce now a practical 
notation that makes it shorter and thus easier to manipulate in convoluted compu- 
tations. For a finite element function, v € V" (or more generally for any Lipschitz 
continuous function v that is C^(int(Jsr)), for each K e denote by A^xv the reg- 
ular part of the distribution — div (AVu), which is defined as a piecewise continuous 
function such that 

(2.3) {A,xv,4>)^ (- div (AVw) , 0) , V,/> e h1(17). 

The operator A^i is sometime refered to as the elementwise elliptic operator, as it 
is the result of the application of — div (AV-) only on the interior of each element 
K £ £^n- This observation justifies our subscript in the notation. We shall write 
the representation ()2.ip in the shorter form 

(2.4) a(z;,0) = (Aeit',</.) + (JH,(/.)s„, V(/. G Hi(f7). 

Let us now recall some more basic definitions that wc will be using. The discrete 
elliptic operator associated with the bilinear form a and the finite element space V 
is the operator A" : Hj(f7) ^ V" defined by 

(2.5) (A%,0„) =a(^^>n), V0„gV", 

for V G Hj(ri). The L2 -projection operator is defined as the operator Pq : L2(ri) — *■ 
V" such that 

(2.6) {P^v, cl>n) = {v, 0„) , V0„ G V" 

for V G L2(r2); and the elliptic projection operator P" : HQ(r2) — > V" is defined by 

(2.7) a(P>,0„) =a(f,0„), V0„ G V". 

The elliptic reconstruction, which we define next, is a partial right inverse of 
the elliptic projection |MN03| . (Notice that a similar operator has been introduced 
by Garcia- Archilla & Titi, albeit with different applications in mind than ours 
IGATOOp 
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2.2. Definition (elliptic reconstruction) We define the elliptic reconstruction oper- 
ator associated with the bilinear form a and the finite element space V" to be the 
unique operator : Hj(fi) Ho(^) such that 

(2.8) a(^"i;,0) = (A"i;,0), V0 e Hi(f}), 

for a given v G HQ(ri). The function is referred to as the elliptic reconstruction 
of V. A crucial property of is that v — S^^v is orthogonal, with respect to a, to 

yn. 

(2.9) a {v - ^"w, (^„) ^ 0, V(/)„ G V". 

From this property and by now standard techniques in a posteriori error estimates 
for elliptic problems |AO00i IBraOli IVer96| , it is possible to obtain the following 
result. 

2.3. Lemma (elliptic reconstruction error estimates) For any v G V" the following 
estimates hold true 

(2.10) - i^li < — \\{A,iv - A"v)hJ + ^ J[v]hl/^ 

a a 

(2.11) ||,^% - v\\ < Ce.2 \\{A,iV - A''v)hl\\ + Cw.2 \j[v]hT 

where the constants Ck.j are defined in [LM06|, Appendix B]. 
We will use the following shorthand 

(2.12) LU^^U. 

2.4. Definition (discrete time extensions and derivatives) Given any discrete func- 
tion of time — that is, a sequence of values associated with each time node tn — e.g., 
([/"), we associate to it the continuous function of time defined by the Lipschitz 
continuous piecewise linear interpolation, e.g., 

(2.13) U{t) Z„_i(t)C/"-i + ln{t)U'"\ for t G /„ and n G [1 : A^]; 
where the functions / are the hat functions defined by 

(2.14) ln{t) (t) - , ,(i), for t G [0,r] and n G [0 : A^], 

Tn Tn+1 

Ix denoting the characteristic function of the set X. The time- dependent elliptic 
reconstruction of U is the function 

(2.15) ujit) ;„_i(t).?^"-iC/"-i + Z„(t).^"f7", for t G /„ and n G [1 : A^]. 

Notice that a; is a Lipschitz continuous function of time. 

We introduce the following definitions whose purpose is to make notation more 
compact: 

(a) Discrete (backward) time derivative 

jjn 

(2.16) aC/" := 

Notice that dV'' ~ dtU{t), for all t G /„, hence we can think of as being 
the value of a discrete function at tn- We thus define dU as the piecewise linear 
extension of {dU"')n, as we did with U. 

(b) Discrete (centered) second time derivative 
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(c) Averaged (L2-projected) discrete time derivative 

rrn pnrrn — 1 

(2.18) dU"" := P^dU" ^ 2 ^ VnG[l:iV]. 

Tn 

The reason we introduce this notation for is that dU^\ in general, does not belong 
to the current finite element space, V", whereas dU" does. 

2.5. Remark (pointwise form) The discrete elliptic operators A" can be employed 
to write the fully discrete scheme (|1.15p in the following pointwise form 

(2.19) dU'^ix) + A"C/"(a;) = P^rix), Va; e il. 

Indeed, in view of 9?/" + - G V", (fTTSll . and jl^l), we have 

(2.20) 

<Jy4"[/" + dU" - P^'f " , = + ac/" - Po'7", P^V^ 

for any (f) £ Hj(f}). Therefore the function dU"^ + A^U"^ - P(5V" vanishes. 

2.6. Error equation. Let us consider the (full) error, the elliptic reconstruction 
error and the parabolic error which are defined, respectively as follows 

(2.21) e^U-u, 

(2.22) e^uj-U, 

(2.23) p = uj-u. 
We have the following decomposition of the error 

(2.24) e = p-e. 

We can also readily derive the following error relation for the parabolic error in 
terms of the reconstruction error and the reconstruction itself jLM06| : 

{dtp{t),cP) + a (pit), 4>) = {dte{t), 0) + a (c.(i) - c.", </.) 
^^■^^^ + (P„"C/"-i - C/"-\ 0) + (Po"r - fit), 0) 

for aU G Hj(f7), t e /„ and n e [1 : N]. 

3. The duality approach to a posteriori error estimates 

In this section we expose the combination of the elliptic reconstruction and the 
duality techniques in the simple situation where the scheme is discrete only in 
space. This should help understand the main ideas that we will employ for the 
fully discrete case in 21 

3.1. Spatially semi-discrete scheme and continuous time notation. In this 
section, since we will be dealing with the space semidiscrete scheme only, we will 
use the same symbols introduced for the fully discrete scheme in fj2l albeit in their 
semidiscrete analog by dropping the index n. The notation now introduced is valid 
only mOand it is important not to confuse time-dependent functions, such as U, 
to, e, e and p, to be introduced next with their fully-discrete analogs introduced 
in [21 Specifically, let V be a given (time-invariant) finite element space, we define 
?7 : [0, T] — *■ V to be the solution of the semi- discrete Galerkin finite element 
equation 

C/(0) := lu{0), and 
^^■^^ {dtU{t), $) + a {U{t), $) = {fit), $) , V$ e V, i e [0, T], 
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where the operator / is a suitable interpolation or projection operator from H^(r2), 
or L2(r2), onto V. Following on this definition we define the error at time t to be 
e{t) := U{t) — u{t) and the semi-discrete elliptic reconstruction to be ijj(i) ■.= MU{t)^ 
were 3£ is the elliptic reconstruction operator associated with the space V (compare 
with 12. 2( ). In analogy with the fully-discrete notation in 12.41 we define the elliptic 
reconstruction error e := lj — U and the parabolic error p := U — u. 

We observe that while in the semi-discrete setting one assumes the discrete space, 
V, to be invariant in time, in the fully discrete setting, analyzed later in|H we will 
take into account the possibility of the discrete space to change, with respect to the 
timestep. For instance, in an adaptive mesh refinement scheme the space change 
derives from the mesh's modification from a time to the next. 

The notation here introduced leads to the semi-discrete analog of the error equa- 
tion ((T^ 

(3.2) {dtp{t)A) + a (pit), 0) - (dteit), 4>) + (PJit) - f{t),<j>) , 

for (f> S Hj(fi),t G (0,T]. This equation consitutcs the starting point for the 
estimate derivation. 

3.2. The dual solution. The concept of dual solution, with that of elliptic recon- 
struction, consitutc the two main analytical tool of this paper. For each s < T, we 
define the dual solution to be the function 

(3.3) z{x,t;s) = Zs(x,t), for x Cz and < t < s, 

which satisfies £ L2(0, T; Hj(rj)), dtZs E L2(0, T; H"^(fJ)), and solves the follow- 
ing backward parabolic dual problem: 

^^^^ -(5tZ«(t),0)+a(0,z.(t)) = O, V0eHj(l]),te [0,s), 

Zs{x, s) = p(x, s), Vx e O 

for each s £ [0,T]. Notice that (j) can be taken to be time dependent, with the 
appropriate differentiability properties. 

The dual solution enjoys stability properties which we will use in the sequel. An 
immediate property is the usual energy identity 

(3.5) ||z,(t)||V2 J|z3|^ = ||p(,s)||^V^e[0,.]. 

A more intricate property is given by the following result. 

3.3. Lemma (Strong stability estimate |EJ91i Lcm. 4.2]) Suppose that ft is a 
convex domain. Then for each s E [0, T] it holds 

(3.6) 

\\dMt)f {s t) dt, Ih div {AVzs) {t)f is - t) dt^ < i ||p(,s)||^ . 

3.4. Error analysis in the spatially semi-discrete case. Integrating by parts 
in time, problem p.4p implies 

{p{s),^{s))^{zs{s),<l>{s)) 
^^■^^ = {zs{0), m) + I {dtm.zs{t)) + a m),zsit)) dt 



for aU (f> e L2(0,r;Hj(r2)) such that dt(t> e L2(0, T; H"^(J1)). 

Taking (p = p, using (j3.2p and assuming Pq/ — f = momentarily — in the proof 
of Theorem 14.21 we remove this assumption — we obtain 



(3.8) \\p{s)f = (p(0), z,(0)) + / {dte{t), z,(t)) dt. 
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The first term on the right-hand side, is easily estimated, with Lemma 13.31 in mind, 
as follows 

(3.9) (p(0),z.(0)) < ||p(0)||sup||z,|l. 

[0,.s] 

As for the second term on the right-hand side of (|3.8|1 we have the choice of two 
different ways for estimating it. 

(a) A direct estimate yields 

(3.10) / (9te,z,) <sup||z,|| / \\dte\\ . 

Jo [0,s] JO 

Notice that the term dt^ can be estimated via elliptic a posteriori error estimates 
because it is the difference between dtU and its reconstruction MdtU ~ dt^U. 

(b) A less direct estimate, that would avoid the appearance of time derivatives 
in the estimator, is obtained by integrating by parts in time first 

(3.11) / {dte,zs) = (e(s),2,s(s)) - (e(0),2.(0)) - / {e{t),dtz,{t)) dt. 
Jo Jo 

The last integral can be then bounded as follows 

,it)dtz,{t)dt< r ^^^\\dtz,{t)\\v^tdt 

Jo Vs-t 

Unfortunately this bound may not be useful. Indeed, for the first integral in the 
right-hand side to be finite it is necessary that e{t) = o(l) at i = s. This means 
that the error between the discrete solution and its reconstruction should at least 
vanish at s. Heuristically this can be interpreted as the mesh having to become 
infinitely fine as time gets closer to s: an unrealistic option. 

To circumvent this difficulty, without totally sacrificing j|e|| to ||9fe||, we com- 
promise between approach (a) and |(b)| by following through from ()3.8p as follows: 
fix r G (0, s) (think of it as a close point to s), split the integral and integrate by 
parts in time 

\\pis)\f^{zM,Pm + (r+ n idte,z,) 



(3.13) 



(z,(0),p(0)-e(0)) + (z,(r),e(r)) - / {e,dtZs) + {dte,z,) 



<sup|jz,||(|je(0)|l + ||e(r)||+ / 

[0,s] 



\\dtzs{t)f{s-t)dtj \ I ^^dt 

The stability estimates p.Sp and (j3.6p imply that 



1/2 / rr llg(i)||2 



1/2 



/ r 2 \ 

(3.14) ||p(.s)||<||e(0)|| + ||6(r)||+ / \\d,e\\ + l{ f dt\ , 







We have essentially proved the following result. 



10 



OMAR LAKKIS AND CHARALAMBOS MAKRIDAKIS 



3.5. Theorem (Semi-discrete duality-reconstruction a posteriori error estimate) 
Suppose that f{t) G V, for t G [0,T], and that there exists an a posteriori cUiptic 
error estimator function •, -l such that 



(3.15) 



\\v^3gv\\ < ^[v,f,L2in)l forve Hj(f}), 



then the error occuring in the semi-discrete scheme p.ip obeys the a posteriori 
bound 



(3.16) 

where 
(3.17) 



sup \\U{t) -u{t)\\ < \\U{0) - u{0)\\ + L{s,r)su-pr?[U,W,U{n)] 
te[o,s] lo,r] 

+ {s-r) snp>^[dtU,Y,L2{n)] 

[r,s] 



L(s,r) 



2 V s — r 



Proof Fix r and s, use (|3.14p with s' < s instead of s, r' = s' — {s — r) instead of 



3.6. Corollary (Semi-discrete duality-residual a posteriori estimates) If is a con- 
vex domain in M'' and f (t) G L2 (il) for each i e [0, T] , then the foUowing a posteriori 
error estimate holds 

sup|lJ7-Mi| < ||?7(0) - u(0)|l 

+ L{s,r) sup (C3 Wh'^iAU - AhU)\\ 

[0,r] 



(3.18) 



C5 



h^/^J[U] 



+ (s - r) sup {C3 \\h^AdtU - AhdtU)\ 



h{Pj-f) 



Proof Sketch: apply the Theorem and add the / part. 



4. Estimates for the fully discrete scheme 

Having understood the main idea of the techinque in the last section, we turn 
our attention to the analysis of the fully discrete scheme (|1.15p . For convenience, 
we switch notation slightly and use the symbol U , even without the superscript n, 
for the fully discrete solution and its pieccwice linear interpolation now. 

We introduce first some extra notation to be used in this section. 

4.1. Definition (Error indicators and estimators) Suppose an a posteriori elliptic 
error estimator function is available, introduce the following (time- local) 

S-hased spatial error estimators 

(4.1) £„ =.<f[C/",V",L2(f^)], 

(4.2) 77„ = <f[9tt/",V"nV"-\L2(r!)], 

and the time estimator 
(4.3) 

On = \\A"-^U''-^ - A^'U' 



Pop - dU^ - A^U" 



for n = 1, 
for ne[2: N] 
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and the logarithmic time accumulation coefficients 



(4.4) 

(4.5) a„ 
where 

(4.6) A(x) 



^rilog(^^^) forne [1 :iV-l], 
\ i for N; 

^ (t^) - ^ {-T^) ' for n e [0 : 71 - 2] , 
A ( ) - 1, for n = iV - 1, 



(1 + 1/x) log(l + x) for \x\ e (0, 1), 
1 for .T = 0, 



which is an increasing function of x. We observe that the functions X{x) — 1, 
1 — A(— x) and A(a;) — A(— y) are positive for {x, y) E (0, 1)^, a fact that makes the 
coefficients a„ to be positive. 

FinaUy we introduce the data approximation estimators 

(4.7) f3n= r Wr-fm di; 

and we notice that for / regular enough we can replace /?„ by the right hand side 
of the inequality 

(4.8) r \\r - fm dt < rn. 

4.2. Theorem (General duality a posteriori parabolic error estimate) Let u he the 
exact solution of {232), (^")ne[i Ar] corresponding (fully) discrete solution of 
il.l5]) and (oj") = (.^"C/) the elliptic reconstruction ofU, as defined by then 
the following a posteriori error estimate holds 

(4.9) 

CN~1 \ 
Tl = / 

/TV \ AT 

V 2 



L2(0)) ' 
n=l 



y^||7Jv/iw|| + ^E ^" j 



1/2 



wiere tie various error estimators and the coefficients Sn, Vn, o-n, bn, Pn, 3.nd 7„ 
are defined in ^4.1^ 

An alternative time error estimator is given by replacing the fourth summand 
on the right hand side by i4.9\) by 



N \ 1/2 

2 



(4.10) E^"^^' 
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4.3. Corollary (Duality a posteriori full error estimates) With the same notation 
as in \4.2\ we have 



N 

K - u{T)\\ < - u{0)\\ + J^Wj^h^W + ^ r„ |19t/|lL,(z„ 



;L2(n)) 



("^•11) +\/l + log — ( max £„ + 2 max hnhll 

" Tat \n6[0:Ar] ne[l:Af-l]" ' 



+ - max (\\U"-^ 

2 nG[l:Ar] ^" 



C/1|+77„)). 



Furthermore, we can chminatc all the terms e„ from the estimate by replacing the 
last term by 

(4.12) max 6'„. 

4.4. Remark (comparison between Theorem 14.21 and Corollary 14. 3p Corollary 14.31 
has a simpler estimate than Theorem 14. 2 1 in that it involyes less terms and does not 
require as much memory. Notice howeyer, that from an error bound view-point, 
the Theorem's tighter bound may be more effective as the time accumulation is 
not as strict as in the Corollary. This is especially true in problems, typical in the 
parabolic setting, where the initial error may be very big and gets damped with 
time. 



5. Proof of the main results 

As with the semi-discrete case of ^ the starting point is the fully discrete analog 
of (|3.8p . which is readily obtained from p.25p and (j3.4|) . and reads 

(5.1) ||p(r)||^ = (p(0), zt(0)) + r (dteit), ZT{t)) dt 

Jo 

N 

/ a{iuit)-Uj\zT{t))+hn,ZTit))+lr~f{t),ZTit)) dt. 

where 

pnTrn—l jjn—1 

(5.2) 7„ := -° + PS~r - r = (^^0^ - I)(r-^C/"-^ + /")• 



5.1. Space error estimate. The first two terms arc estimated, similarly to ()3.13p . 
as follows 

(5.3) (p{0),ZTm+ [ {dte{t),ZT{t)) dt< 
Jo 

T ft 2 \ -'^/'^ 

\\p{m I l|e(0)|| + ||e(t^-i)|| + ^ ||a,e|| + ^ N ""^^^ dt j 
To proceed we observe that 

i-T 



(5.4) / listen ||e^-e^-i II <77Ar. 

Af-1 
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and that, by convexity and linearity, 

2 



(5.5) 



Thus we obtain 
(5.6) 



T-t J a T-t 

Tl = "^^ Tl = 



(p(0),zt(0)) + [ {dte{t),ZT{t)) dt < ||p(r)|| I ||e(0)|| + ( E 1 + 



We notice that the former estimate imphes the more traditional one [EJ91j 

(5.7) ^0 T-t -neiO:N-n' [ 

= (l+4X(T,.._,)^)„max_^,lk"f 

where _L(T, t^v-i) is the logarithmic factor defined in p.l7p . 

5.2. Time error estimate. The third term in (|5.ip . which accounts mainly for 
the time error, can be bounded as follows 

(5.8) 

E / a{uj{t)-oj",ZT{t)) dt<Y^ / ||a)"-i-w"||/„_i(t)||-div(AVzT)(Oll 

<^iip(^)ii (EK-^--ir/"^^Y^dt^ ' 
<^ib(T)ii (iK--^-^f +Eii-"--""'iri°g(^^^)) ' ■ 

Notice that we do not need to deal with the last interval (tAr_i, T] as a special case 
since lN-i{t) = 0(T — t) compensates for the singularity of 1/(T — t). The terms 
||a;"~^ — cij"!! appearing in this estimate still need to be estimated, as there is no 
explicit knowledge of the reconstructed functions uj" = ^"[/". These terms can be 
dealt with in two different ways. 

(a) One way to estimate these terms is given by: 

- uj"\\ < II [/"-I - ;7"|| + ||cj"-i - - u"-^ + C/"|| 
^^■^^ = ||t/""'-f/"||+T„||ate"|l < ||c/"-'-t/"||+7M, 

for all t E In- Thus we obtain the estimate 
(5.10) 

N »f„ / N \ 1/2 



E / aiujit)-u;",ZTit))dt<\\piT)\\lY,bn{\\U"'-'-U''\\+rin]'' 

,1=1 \n=l 
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(b) Another way to estimate ||a;" ^ — a;"!! consists in using again the definition 
of eUiptic reconstruction and the Poincare inequahty as foUows: 

(5.11) = C2,i - - t^") 

< C2,i - A"C/"|| - 

thus obtaining 

(5.12) " ^"|| < C'2,1 - A"C/"|| = 0„. 
Therefore 

(5.13) 5^ / a (c.(0 - a;", ^^(t)) d< < ||p(T)|| ^ fo^^f, 

Also here, the estimator simphfies if we relax the bound by taking the maximum 
norm in time 

(5.14) 

V /*" a{uj{t)-iu",ZTit)) dt<\\p{T)\\ max - |L / f i + L(r, tw_i)2 

5.3. Remark The above bound could be derived differently by using not the strong 
stability, but rather the definition of elliptic reconstruction directly. We refer to 
our companion paper for a detailed estimate of this type [LM06| § 3.5]. 

5.4. Data approximation and coarsening error estimates. We bound now 
the third term in (15.11): 



^ rU^ N ,t. 



^-1 /•*„ rT 



^Y. r h^^lW \^T{t)\^ dt+ I hNhNW \zT{t)l dt 

(5-15) , . \ 1/2- 

<\\\p{T)\\ (^\\^MhM\\+(Y.\\lnhl\\^\0g(^'^'^-' 



\piT)\ 



1/2N 



Here we have used the fact that is convex in order to apply the estimate 
(5.16) |zT(t)|2 < ||-div(AVzT)(t)||, 

and then apply the strong stability estimate (j3.6p . As with the space and time 
estimates, this estimate can be simplified, with some loss of sharpness, as follows 



(5.17) f" ^ln.ZT{t)) dt 

n=l-^*"-i 



< \\piT)\\ \ -z-\\lNhN\\ + L{T,tN-i) max 7„/i,, 

\ V 2 ne[l:N-l] " 

The fourth term in ()5.ip can be bounded in two different ways depending on 
which definition for /" is taken. 
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(5.18) 



(a) If /" = /" then we can proceed as follows 

n=l-^*"-i 



E / (/" " /(^)' ""Tit)) dt < E max IIztII / /" - f{t) 



At 



N 



< 



\piT)\\Y.f^- 



(b) If instead we have /" = J^" ^ /(t)di/T„ then we can exploit the orthogo- 



nality in time and write, for each n [I : N — 1] 



r-f{t),ZT{t)) dt 



(5.19) 

By noticing that 
(5.20) 



max||zT(i) — ZT{tn-i)\\ ~ max 



r^f{t),ZT{t)~ZTitn-l)) dt 



<max||zT(t)-^T(t„-i)|| / r-fit) 



dt. 



/ dsZT{s)ds 



< max / \\dsZTis)\\ ds < / \\dtZT\\ 



< log 
= 2bl/' 



T — tri-l 

T~t„, 



1/2 / „t 



1/2 



\\dtZT{t)r {T - 1) dt 

\ 1/2 



\\dtZTit)\\^{T-t)dt 



t„-i 



Summing up and using the strong stability estimate (j3.6p we obtain 



(5.21) / (/" - /W' ^^(*)) < /3a^ + 2 E ^"A 



W-l N 1/2N 

2 



Like earlier estimates, this estimate can be further simplified, by taking the maxi- 
mum and slightly relaxing it, into 



(5.22) E /" {r-f{t).ZT{t))dt 

n=l 



< |1p(T)|1 Pn + 2L(T,iAr_i) max /?„ 

ne[l:7V-l] 



Replacing the estimates ((S^ . ((CTI)) / ((O^ , ((OK)) and ^^A^ into the relation 
()5.ip we obtain the result of Theorem [ 



5.5. Proof of Corollary 14.31 The splitting p = e + e and some simple Holder 
inequality manipulations on the summations on the logarithmic coefficients yield 
the result. 
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6. Sample application: residual estimates in space 



6.1. Definition (residuals) The residuals constitute tlie building blocks of the a pos- 
teriori estimators used in this paper. We associate with equations (|1.12|1 and (|2.19p 
two residual functions: the inner residual is defined as 

i?° := Ael[/° - 

(6.1) _ 

i?" := AeiL/" - .4"[/" = AeiJJ" - Po"/" + 5C^", for n£[\:N], 

and the jump residual which is defined as 

(6.2) r := J[C/"] = lVt/"l . 

With definition fj2Tl]in mind, the inner residual terms can be written explicitly as 



(6.3) (i?",0>= (-div(AViO-Po"r 



Ijn _ pnjjn-1 



K 



We can now introduce, for n G [0 : A^], the elliptic reconstruction error estimators 



(6.4) en:=C6,2||/i„i?"||+Cio,2 
and, for n £ [1 : N]^ the space error estimator 



(6.5) 



Vn 



6,2 



C 



10,2 



c 



14,2 



6.2. Conclusion. With these definitions. Corollary 14.31 includes the estimates of 
Eriksson & Johnson |EJ91| . but with no unrealistic assumptions on the meshes 
from time to time. 

Note also that a direct application of Theorem 14.21 provides finer estimates with 
respect to time accumulation; these are especially helpful in situations where the 
error (on a time-invariant mesh) decreases with time. 
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